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Abstract. We develop an approach to the theory nonholonomic relativistic stochastic processes in curved 
spaces. The ltd and Stratonovich calculus are formulated for spaces with conventional horizontal (holo- 
nomic) and vertical (nonholonomic) splitting defined by nonlinear connection structures. Geometric models 
of relativistic diffusion theory are elaborated for nonholonomic (pseudo) Riemannian manifolds and phase 
velocity spaces. Applying the anholonomic deformation method, the field equations in Einstein gravity 
and various modifications are formally integrated in general forms, with generic off-diagonal metrics de- 
pending on some classes of generating and integration functions. Choosing random generating functions 
we can construct various classes of stochastic Einstein manifolds. We show how stochastic gravitational 
O | interactions with mixed holonomic/ nonholonomic and random variables can be modelled in explicit form 

and study their main geometric and stochastic properties. Finally, there are analyzed the conditions when 
1 non-random classical gravitational processes transform into stochastic ones and inversely. 

a' 

PACS. 02.50.Ey Stochastic processes - 02.90.+p Other topics in mathematical methods in physics 
02.40.-k Geometry, differential geometry, and topology - 04.20.Jb Exact solutions 

m ■ 

> 1 1 Introduction 

In recent years stochastic methods and the theory of (relativistic) diffusion on Lorentz manifolds, and (for instance) 
Finsler/ supersymmetric generalizations, re-attracted considerable interest and applications in various directions of 
mathematics, gravity and modern cosmology and astrophysics. During the last decade a lot of papers have been devoted 
to such issues (there are many ways and approaches); for reviews and original results in physics and mathematics, see 
[2PP||31|rjl|5pi|II^] and references therein. 

In a series of our works [T7lll8ll2T] , we proved that the gravitational field equations in various theories of gravity 
(general relativity and extra dimension models, noncommutativc, nonsymmetric, Finsler etc generalizations) can be 
integrated in very general forms using the so-called anholonomic deformation/ frame method. The idea of the method 
is to use some distortion tensors completely defined by a given, in general, off-diagonal (which can not be diagonal- 
ized by coordinate transforms) metric structure and uniquely deforming the Levi-Civita connection into "auxiliary" 
1 linear connections which are also compatible and uniquely defined by the same metric. For such a connection and 
with respect to some classes of nonholonomic frames of reference, the system of Einstein equations decouple split 
into conventional sub-systems which can be integrated "step by step" in very general forms. Imposing additional 
nonholonomic constraints on integral varieties, we can generate exact solutions in the general relativity theory. 

In explicit form, the type of exact solutions depend on the class of generating and integration functions and 
considered boundary /initial conditions, on prescribed smooth classes of coefficients and sources of equations and on 
existing spacetime local and global symmetries for a chosen topology. We have to involve additional physical arguments 
in order to select a "physically important" exact solution with deterministic or random properties. The bulk of exact 
solutions in gravity theories were constructed for smooth and/or singular classes of metrics and there were analyzed 
possible perturbative deformations by quasi-classical quantum contributions and fluctuations of classical and quantum 
matter sources. 

In this work we study new classes of stochastic gravitational processes modelled as exact solutions of the Einstein 
equations (of vacuum type, with cosmological constant and/or with various types of stochastic sources). There are 
developed certain methods for the theory of stochastic calculus and relativistic diffusion on curved spaces following 
the formalism of moving frames with associated nonlinear connection structure and respective constructions from 
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the geometry of nonholonomic manifolds. We prove that such new classes of stochastic solutions in gravity are de- 
rived by integrating in general form certain systems of PDE with separation of equations^ and stochastic generating 
functions and by imposing corresponding nonholonomic (equivalcntly, anholonomic, or non-integrable) constraints on 
gravitational vacuum and non-vacuum interactions. Our aim is also to formulate special and general relativistic 
nonholonomic diffusion models and apply such stochastic techniques in modern gravity theory. 

With respect to exact sure and stochastic solutions of the Einstein equation, we have to address to the general 
problem of transition to stochastic phase of solutions for a system of nonlinear partial differential equations 
(PDE) on curved manifold/ bundle spaces is the following: 

— If and under what conditions can an existing "non-random" solution (depending on certain smooth classes of 
generating and integration functions) be further transformed into a stochastic solution of the same system of 
equations? 

— What types of nonholonomic constraints, boundary /inicial conditions, and corresponding data for generating/ in- 
tegration functions will transform a nonlinear dynamical/ evolution classical, or quantum process, into a stochastic 
one, and inversely? 

In brief, this can be equivalently posed as a problem for formulating some criteria for generating nonlinearly 
stochastic processes (for certain well defined conditions, diffusion). In general, the solutions of systems of PDE can be 
with mixed sure and random variables. 

We shall construct exact stochastic solutions of Einstein equations in explicit form for vacuum and non-vacuum 
gravitational configurations. Such spacetimes are described by generic off-diagonal metrics, the nonlinear gravita- 
tional interactions being subjected to certain types of nonholonomic (equivalently, anholonomic, or non-integrable) 
constraints/ conditions which may give rise, or not, to random gravitational processes. A cosmological constant, with 
possible local anisotropic polarizations, and matter sources present additional issues to be addressed^ 

A partner work [T| contains developments of the results of this paper for stochastic Ricci-Finsler flows and gravi- 
tational interactions/ evolutions modeling porous media and self-organized criticality. 



1 it should be not confused with the method of separation of variables 

2 For non-experts on probability and stochastic calculus (i.e. researches skilled in differential geometry and applications), we 
remember some main concepts and ideas which are important for definition of relativistic diffusion theory (see details in [221 
I23II24||2"5] ): Physicists have abandoned determinism as a fundamental description of reality and began developing probabilistic 
(that is, stochastic) models of natural phenomena long before quantum mechanics and physics. Classical uncertainty preceded 
quantum uncertainty (in our approach derived from solutions of nonlinear systems of equations). P. Langevin, was the first 
to apply Newton's second law to a "Brownian particle" on which the total force included a random components. The time 
evolution of a random variable is called a random or stochastic process. Thus X(t) denotes a stochastic process. The time r 
evolution of a sure variable is called a deterministic process and denoted x(t). 

Random and Sure Variables: A quantity that, under given conditions, can assume different values is a random variable. It 
matters not whether the random variation is intrinsic and unavoidable (for instance, a consequence of quantum effects) or an 
artifact of our ignorance. A random variable is conceptually distinct from a "certain" or "sure", i.e. deterministic, variable. 
The expected value of a random variable X is a function that turns the probabilities P(x) into a sure variable called the mean 
of X. The mean (P) parametrizes the random variable X, but also do all the moments (X n ) (n — 0, 1, ...) and moments about 
the mean ((X — (X)) n ). A continuous random variable X(t) is completely defined by its probability density p{x). 
A memoryless process is a Markov process. The Wiener (Browninan) process, defined by the Markov propagator with a 
"simple" parameter p 2 is the simplest of all Markov processes. The Wiener process is, on its domain, everywhere continuous but 
nowhere smooth. This special property makes the Wiener process dynamical equation a different kind of mathematical object 
- a stochastic differential equation (SDE). It is characterized by a diffusion constant. For a stochastic process X(t), the mean 
function is p,(t) = E[X(t)] and the covariance function is k(t,t') = E[(X(t) — u(t)) (X(t') — fi(t'))]. Gaussian processes are 
stochastic processes defined by their mean and covariance functions. 

A diffusion equation is mathematically equivalent to a stochastic dynamical equation. In gravity theory, random variables 
are introduced when some components of metric are changed by stochastic forces (generating functions) driven by a Wiener 
process. Markov process can be considered but with "some" memory on hyperbolicity (finite speed of light), Lorentz transforms 
in special relativity, and Mach principle in general relativity, as well on Einstein equations etc. This is encoded in the definition 
of Stratonovich integral on curved/nonholonomic spaces. 

3 A nonholonomic manifold is defined by a pair (V,Af) , where V is a manifold and AT is a non-integrable distribution. We 
have to involve in our research certain methods from the geometry of nonholonomic distributions with associated nonlinear 
connection structure (geometry of nonholonomic manifolds) and introduce on Minkowski and (pseudo) Riemannian spaces some 
types of Lagrange-Finsler parametrizations of metrics and connections because only using such variables we can construct 
general exact off-diagonal solutions in gravity and formulate an unified mathematical approach to classical field interactions 
and random processes with mixed holonomic/ nonholonomic / stochastic components. 
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1.1 Related directions 

In order to avoid possible ambiguities with the formalism and terminology used in different approaches, we mention 
here some alternative methods with stochastic/ diffusion in curved spaces and random variables in gravity: 

1. The theory of stochastic semiclassical gravity (in brief, stochastic gravity) is based on the Einstein-Langeven 
equation with additional sources due to the noise kernel. It was naturally constructed from semiclassical gravity 
and quantum field theory in curved spacetimes and non-equilibrium statistical mechanics, see a comprehensive 
review [35] and references therein. For short, stochastic gravity includes also its fluctuations, for instance, of 
quantum fields in curved spaces and explore how the metric fluctuations are induced and seed the structure of the 
universe, how such processes affect the back reaction of Hawking radiance in black hole dynamics and the black 
hole horizons, trans-Planckian physics etc. 

2. The theory of stochastic (diffusion) equations on curved manifolds with local Euclidian signature^ is a well de- 
veloped topic (see e.g. 23,24,25 ), see also related applications to locally anisotropic kinetic processes, gravity 
and astrophysics 11,12 and an alternative approach to anisotropic diffusion [32] . For such constructions, random 
processes were considered on classical (manifolds) inducing possible matter field sources for gravitational filed 
equations. 

Here we also note that various stochastic methods were also applied and developed for quantizing gauge and 
gravitational fields [3~3 l34|[3"5] . see also a recent attempt for stochastic quantization of Hofava-Lifshitz gravity |36) . 
In this work, we shall restrict our considerations only to nonlinearly defined classical (non-quantum) gravitational 
processes generating exact solutions of the Einstein (not obligatory in the form of Einstein-Langeven) equations. 

1.2 Goals of the paper: 

Conventionally, there are three main purposes: 

1. To provide an introduction into the theory of stochastic processes on nonolonomic Euclidian and Riemannian 
manifolds. 

2. To generalize the theory of stochastic equations in relativistic form for nonholonomic (pseudo) Riemannian mani- 
folds in a form allowing to study mixed sure (non-random) and stochastic gravitational processes. 

3. To prove that prescribing stochastic generating and integration functions for "formal" general solutions of Einstein 
equations we define three types (so-called, horizontal, vertical and horizontal-vertical ones) of stochastic (non- 
holonomic) spacetimes. It will be shown how to construct general stochastic gravitational solutions with Killing 
symmetries and extensions to non-Killing ones and formulated the criteria when the non-stochastic solutions may 
transform into some stochastic ones, and inversely. 

We organize the exposition as follows: In section [5] we summarize the Ito and Stratonovich stochastic calculus and 
diffusion theory on nonholonomic Euclidian and Riemannian manifolds. We use a synthesis of geometric and stochastic 
methods originally elaborated for the diffusion theory on nonholonomic vector bundles spaces and in Lagrange-Finsler 
geometry. Section [3] is devoted to the theory of nonholonomic special and general relativistic diffusion. A formalism of 
adapting the constructions to the nonlinear connection structure play a key role in distinguishing nonlinear gravita- 
tional and random gravitational processes and keeping certain similarity with former constructions for phase velocity 
spaces and moving frame method. In section 01 we generalize the anholonomic deformation method of constructing 
exact solutions in gravity in order to include into the scheme the possibility to generate stochastic metrics and non- 
linear and linear connections. We also analyze explicit conditions/criteria when such stochastic gravitational fields are 
of vacuum type, with gravitational polarizations and/or cosmological constants and state the possibility to extend the 
constructions for nontrivial matter sources. A summary of results and concluding remarks are given in section [5l 



2 Stochastic Processes and Nonholonomic Manifolds 

In this section, we provide an introduction into the theory of stochastic differential equations on nonholonomic 
(pseudo) Riemannian manifolds [17]. We assume that the reader is familiar with the concepts and basic results on 
stochastic calculus, Brwonian motion and stochastic processes "rolled" on curved (usually, Riemannian) spaces [221 

mm- 

4 in this paper we shall use equivalently the terms "space" and/or "spacetimes" for certain (pseudo/semi) Riemannian ge- 
ometries, endowed with metrics of necessary signature, nonholonomic distributions with associated anholonomic frames, and 
(auxiliary) generalized connections 
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2.1 Geometric preliminaries 

We consider a 4-d manifold (space, or spacetime) V endowed with a metric 

Z = g afj {if)du a ®du fi (1) 

of necessary signature (±,±,±,±) when local coordinates are parametrized in the form u a = (x t ,y a ), where x 1 = 
(x 1 , x 2 ) and y a — (y 3 — v, y A — y) Indices i, j, k, ... = 1,2 and a, b, c, ... = 3, 4 are for a conventional (2 + 2)-splitting 
of dimension when the general (small Greek) abstract/coordinate indices when a, (3, . . . run values 1, 2, 3, 4. 
For our purposes, we split a metric ([T]) in the form 



-- gij dx l ® da? + h ab (dy a + N£dx k )®{dy b + N b k dx k ) : (2) 



and parametrize as 



r, g = rn(x k ,v) °g l (x k ,v)dx i <g> dx i +i] a (x k ,v) °h a (x k , v)e a ®e a , (3) 
= dv + r] 3 (x k ,v) °w i {x k 1 v)dx\ e 4 = dy A + r)f(x k ,v) °n i (x k 1 v)dx\ 



e 3 



In we consider that gtj — diag\gi — rji °gi] and h ab = diag[h a = rj a °h a ] and iV| — Wi — rjf °Wi and 
Nt = Tii = i]f °ni. The gravitational 'polarizations' r\ a and i]f determine nonholonomic deformations of metrics, 
°g=[ °gi, °h a , °N£] — > ? 'g=[ gi,h a ,N£] and can be defined by functions of necessary smooth class and/or any 
random (stochastic) variables|f| 

We use a boldface symbol V for nonholonomic manifolds (equivalently, spaces) when a Whitney splitting of tangent 
space TV is defined, 

TV = hV®vV, (4) 

with conventional horizontal (h) and vertical (v) subspaces, respectively, hV and uV. Such a geometric object defines 
a nonlinear connection (N-connection) structure N, stated locally by a set of coefficients {N£} with respect to a 
corresponding coordinate basis. It should be emphasized here that we can introduce on (pseudo) Riemannian manifolds 
various types of nonholonomic distributions and N-connections, for instance, with 2 + 2 splitting, following the principle 
of covariance allowing us to consider any equivalently any system of reference/coordinates. 

A N-connection structure N = {N£} is present in @ via iV-adapted frame, e Q , and dual frame, e^, structures 
(i.e. N-elongated partial derivatives, respectively, differentials) 

= di -N?d a ,e b = d b = -^j, (5) 
e- 3 = (e l = dx\ e a = dy a + N?dx l ) . (6) 



Such local bases satisfy nonholonomic relations of type 



[e a , ep] = e a ep - epe a = w 7 Q;3 (u) e 7 , 



with nontrivial anholonomy coefficients w^ 7 (u) , 

which is also related to the concept of nonholonomic manifold jl7U18) . i.e. a manifold endowed with a non-integrable 
distribution (in particular, with N-connection structure). 

On a spacetime V, there is an infinite number of metric compatible linear connections D, satisfying the conditions 
-Dg = 0, and completely defined by a metric g @. A subclass of such linear connections can be adapted to a chosen 
N-connection structure N, when the splitting Q is preserved under parallelism, and called distinguished connections 
(in brief, d-connections)0 A general d-connection is denoted D = (hD,vD), being distinguished into, respectively, 



5 In this paper, we shall consider metrics with local signature (+, +, — , +) and/or (+, +, +, — ) in order to follow conventions 
from our previous results on exact solutions with local anisotropy [171118] . 

6 Such transforms, with deformations of the frame, metric, connections and other fundamental geometric structures, are more 
general than those considered for the Cartan's moving frame method when the geometric objects are re-defined equivalently 
with respect to various classes of systems of reference. 

7 On spaces endowed with N-connection structure, there are used the terms distinguished vectors/ forms / tensors etc (in 
brief, d-vectors, d-forms, d-tensors etc) in order to emphasize that the geometric constructions are adapted to the N-connection 
structure, i.e. preserving h- and v-splitting. 
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h- and v-covariant derivatives, hD and vD. We proved [171118) that the Einstein equations in various gravity theories 
decouple and became integrable in very general forms for the so-called canonical d-connection, D. With respect to 
N-adapted bases §5§ and ([5]), D is computed to have the coefficients T~ > aj3 = {L l - k ,L bk , C* c , Cg c ), for hD — {Lj k , L bk } 
and vD = {C} C ,C% C }, with 

L)k = \d lr (etgjr + e 3 9kr - e r g jk ) , L a hk = e b (N%) + hi ac (e k h bc - h dc e b N^ - h db e c N^) , (7) 
C'jc = ^9 lk e c gjk, C£ c = hi ad (e c h M + e c h cd - e d h bc ) . 

The d-connection D and its torsion f = {T 7 Q0 = f 7 Q/3 - f 7 ^ Q ; f l jk , f ) a , f ^ f a u , f a bc }, where the nontrivial coeffi- 
cients ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ 

rpi ri ji rpi s~ti rpa /~)a rpc tc { A7C\ rpa r~ia r~ia (Q\ 

1 jk — ^jk ^kji 1 ja — 1 ji ~ ~~ 11 ji> 1 aj ~ ^ aj ~ a \ j / ' be — °bc — °c&i K°) 

are completely defined by the coefficients of metric g © following the conditions that Dg =0 and the "pure" horizontal 
and vertical torsion coefficients are zero, i. e. T l - k — and T a bc = 0. 

Any geometric construction for the canonical d-connection D can be re-defined equivalently into a similar one 
with the Levi-Civita connection following formula 

rlp = K P + zl p , (9) 

where the distortion tensor Z 1 a p is constructed in a unique forrrjfl from the coefficients of a metric g a /3- 

2.2 h— and u-adapted Euclidian diffusion 

A nonholonomic manifold with conventional h- and v-splitting is with local fibered structure and similar to a 
vector/tangent bundle enabled with N-connection structure. In this work, we develop the approach in relativistic form 
(for holonomic spaces with local pseudo-Euclidean signature when the curved spacetime V has a tangent space with 
splitting of dimension n + to, for n, to > 2. 

A distinguished Wiener process (in brief, Wiener d-process) of dimension n + to is defined locally by a couple of 
elementary (Wiener) h- and v-processes W"(r) = (W(t), W q (t)) , where r (in particular, we can take r — t to be a 
time like parameter) We consider a random (stochastic) curve on V lifted to the horizontal curve on the frame of 
orthonormalized bundles O(V) related by frame transforms, e a > — e— a ,(u)da — e a a ,(u)e a , to N-adapted bases e Q (JS]) 
when by da — d/du— — (pi — d/dx l ,da = d/dy-) we denote a local coordinate base (if it will be necessary, we shall 
underline indices for coordinate bases; primed indices will be used for coordinates with respect to orthonormalized 
frames of reference). 

2.2.1 Ito d-calculus 

A diffusion distinguished process (d-process) in an Euclidian space is described by a couple of horizontal and 
vertical stochastic differential equations, 

dU a = a% (r, U)W a ' + b a (r, U)dr (10) 

where U = (hU, vU) € R" +m is a stochastic d-process with U(0) = u, for u = {u 13 — (x^ , y c )}, with parameter (time 
like variable, or temperature, r > 0). The given values a", and b a are called respectively the diffusion coefficients 

8 we can see this from explicit formulas 

7" „ l.ab 1 r) a 7 i ^ O c h rJ i ~ ih C<3 7 a + B ab T^ c 7 i — ^ Cl a h ni* J_ ~ ih C'i V i — fl 

^jk — ~^jbyik n — — Hjki ^bk — ^ 3 k cb " ~ "f* U hb> ^bk — —cd 1 kbi ^kb — s 'jk n cby +~jk ^hb' Z 'jk — U J 

ij r n 

rya _ — > — >ad rpC rya p. ryi ij rpC j . rpC i 

^jb — — — cb J-jdl ^bc — U, Zj a b — — yljallcb + -tjb'icaj , 

for E% = - g jk g ih ) and ± Sf d = \{8 a J b d + h cd h ab ) 

9 In Euclidian space, the coefficients of a (n + m)-dimensional such Wiener process dW a = W (t + At) — W a (r) are defined 
for the probability density V(W a ) = exp (- ^X-^ ) when the ex P ectations O^™) = and (W a { 1 t)W b ( V+ 2 r)) = 
p{ 2 t) 8 al3 , for S al3 being the Kronecker symbol. 
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and the drift coefficients. It is possible to transform (fTU)) into an integral equation U" = Uq + J cr^, (<;, U)5Wf + 

o 

T 

/ 6 a (<T, U)cfc, i.e. the Ito stochastic integral adapted to h- and v-splitting, defining an Ito process as a Markovian 
o 

process (see details in [23,24,25,8 ). In the above formulas we write 5W" instead of dW" in order to emphasize that 
the approach is with N-elongated partial derivatives and differentials, and ([5]), instead of usual ones. 

If U" is an Ito process, then a function = / (U") is also an Ito process when the Ito N-adapted formula for 
stochastic differential 6f = Af with associated diffusion d-operator A = hA © vA, 

n+m 

A = | E { 2 CT «' (r ' u)ct «' (t ' u) (e ^' + f + <r (r ' u)ct «' (t ' u)e ° e6/ + foQ(r ' u)Gq/} ' 

a' = l 

n 

where, for instance, hA = f X) {h a a' ( r > ^) a a' (r, U)(e i e :) + eje^f + 6 l (r, U)ej/}. Such an operator contains the 

a' = l 

second derivative and the Ito stochastic d-differential Sf does not satisfy usual properties for linear operators which 
makes the theory more sophisticate. 

For any stochastic N-adapted process U" S R" +m , we can introduce the probability density function <p(r, u) and 
determine how it evolves with time/temperature parameter r. This function allows us to compute the probability Pr 
that a realization of a set of variables U\, ...U n + m falls inside a domain IA in the (n + m)-dimensional space of such 
variables, 

Pr(U?eW):= J <j)(T,u)dx 1 ...dx n 5y n+1 ...5y n+m . (12) 
u 

We can calculate the expected value U £[/(U T )] for any function^ / of U, 

u E[f(U T )\ := J f(u)cj>(r,u)Su\..Su n + m . (13) 

If we define /(r, u) := u E[f(XJ T )}, such a function is subjected to the condition (in literature, it is called the 
Focker-Plank, or the forward Kolmogorov, equation) 

a r /(r,u)=A/(r,u), (14) 
/(0,u) = /(r,u), 

where d T = d/dr and the d-operator A is defined by (ITTj) . Physical applications are usually considered following 
exact/approximate solutions of such equations F*l 

Finally, we provide the formula A f = lim T ^o — E ^^ u ^^^ u "> f defined for any suitable function /, where u = U r= o 
is the initial point of N-adapted stochastic process U r . Under the conditions that above integral formulas hold true, 
we can say that a diffusion generator/d-operator A, and its adjoint *A, of U r is associated to an Ito d-process. 



2.2.2 Stratonovich d-calculus 

There is an equivalent reformulation of the stochastic calculus by using the Stratonovich integraQj which is more 
convenient for various curved spacetime generalizations. In N-adapted form, we write 

dU a =5$(t,U) o5W a ' +b a {T,V)dr, 



in physical applications, we can consider a smooth class or any class for which a corresponding integration procedure is 
defined 

11 Similar constructions are possible for the (Hermitian adjoint) d-operator *A of A (the formulas can be proven using (I12|l 
and JT3)), d T (f>(T,u) = *A0(r,u); for arbitrary function /(u), / / (u) * A(/>(r,u)du 1 ...du n+m = J f(u)d T (p(T,u)du 1 ...du n+m , 

u u 

for *A=§ £ {i^,(r,UX(r,U)(e l e J +e J e l )/ + <,(r,UK,(r,U) ea e i) /-e Q b Q (r,U)/}. 

a'=l 

12 even such integrals do not result in Markovian processes 
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where we put " ~ " and " o" in order to not confuse this interpretation with that given by (|10j) for the Ito approach. 
To get equivalent formulations of stochastic calculus is possible if we identify 

n+m 

5=2, (r,U) = <,(r,U), b a (r,V) = 6«( T ,U) - | £ erf, (r, U)e^« (r, U). 

a'=l 

Such a re-definition of drift coefficients results in a linear (on partial derivatives) operator for the associated diffu- 
sion/generator d-operator (compare, for instance, with (jllj) where there are contained second partial derivatives), for 
the Stratonovich interpretation written in the form 

n+m 

A = | J2 L «' L «' + L °' ( 15 ) 

a' = l 

where h a > — & , (r, t^e^ and Lo = 6 q (t, u)e^ are called the fundamental (for a diffusion d-process) d-vector fields. 
The associated N-adapted Focker-Plank equation (compare with p4|l) in the Stratonovich approach is 

n+m 

9 t 0(t,u) = - e /3 crf,(r,u)e 7 [ct^, (t, u)0(r, u)] - e a b a (T, u)0(t, u). 

a' = l 



2.3 Diffusion on nonholonomic manifolds 

The theory of stochastic differential equations on a nonholonomic manifold V enabled with a metric compatible 
d-connection can be constructed similarly to the case of n + m dimensional Ricmannian spaces if we work with N- 
adapted frames and co-frames §5$ and ©■ We use N-adapted variables/coordinates r = (u, e) = (w",e ,3 g) 6 O(V) 
and their frame/coordinate transforms. The infinitesimal motion of a smooth curve it a (r) £ V is lifted naturally to 
that of 7" (r) 6 O(V) using the ordinary differential equations for N-adapted to ((H) parallel transport, 

<5u Q = e a a ,(uP)5j a ' and <Se Q Q ,K) = -f %K)e" a ,K)fa^ 

where we use symbols 8u a , Se a a , etc instead of respective du a , de a a ,..., in order to emphasize that in our constructions 
we use N-elongated partial derivatives and differentials. The coefficients of the canonical d-connection T ^ are 
computed following formulas ([7|). 

A curve r(r) = (u(r),e(r)) is considered as the horizontal lift (it should be not confused with the h-splitting 
considered in previous sections) of a curve u(r) to the nonholonomic bundle O(V) (modelled by atlases of carts 
covered with Euclidean spaces R( rl +' m ) anc [ a "horizontal" curve 7" (r) in the tangent spaces, which can be 

identified locally with an Euclidean spaces R n+m ). Using such lifts, we can define stochastic differential equations 
on nonholonomic manifolds when the fundamental Wiener processes are associated to h- and v-components and 
corresponding Euclidean carts. The corresponding stochastic integrals are defined in the sense of Stratonovich on any 
such open regions of V and O(V) when the canonical realization of multidimensional Wiener processes is used and 
S"f a (r) — > 6W a . The stochastic differential equation describing N-adapted diffusion on a nonholonomic manifold is 

Su a = e a a , (r) o 5W a ' + A a (r)dr (16) 
5eV(T) = -r%(r)eV(r)o^, 

where the components of d-vector A Q (r) = (A l (r), A a (r)) are introduced additionally in order to model various types 
of stochastic processes and take into account possible external forces. Here S a 13 e a a , (u(r))e /3 ^, (u(t)) = g™' 3 with S a 13 

taken as the flat Euclidean metric splitting as S a 13 = (^5 l 3 ' , 5 a b . 

Extending on 0(V) the definition of fundamental d-vector fields from (Tl"5|) . L a , — > °L Q / and Lo — > °Lo, where 

°l a , = e a a ,e a - r%K)e^eV-^-. ° L o = A«(r,U)e a - T^A? e>(r)-^-, 

e p, e p, 

we provide a nonholonomic generalization (horizontal lift of the diffusion d-operator V A) of the diffusion generator 
A on the bundle of orthonormalized N-adapted frames, 

n+m 
a' = l 
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For any projection of a function / in O(V) to V (when, for instance, /(r) = /(u, 0), r = (u a , e^g)), we can write 

°A/(r) = v A/(u), 

where V A = |^ e" Q ,e Q (e^,e p ) + A? e p = P - A + e p (18) 

a.' 

and A = ig Q ' 3 [e a e + e p e a + (f u a0 + f% a ) e „] (19) 

is the canonical Laplace-Beltrami d-operator defined by the inverse coefficients g"^ of d-metric ([2]) , the canonical 
d-connection f" af3 Q) and N-connection N = {N%} (g]). 

The operators V A (|18[) and A (TIT)]) allows us to formulate a generalized Kolmogorov backward equation on a 
nonholonomic manifold V, 

t /(t,u) = V A/(r,u), (20) 
/(0,u) = /(u). 

Since the canonical Laplace-Beltrami d-operator is self-adjoint, A = *A, we can construct a self-adjoint * ( V A^ 

to V A as we explained above in footnote 111! As a result, the corresponding generalized Fokker-Planck equation on 
V is 

d T F = -^=e^J\^\A"F) + ^AF, (21) 

Vl8a/3| V 2 

where F — F(t, 1 u;0, 2 u) is the transition probability with the initial condition F(0, 1 u;0, 2 u) = S( 1 u— 2 u) for 
any two points 1 u, 2 u 6 V and adequate boundary conditions at infinitvP*! Finally we note that we get usual 
evolution/diffusion equations in nonholonomic curved spaces, d T F — f AF, if we impose the condition that the 
divergence in (j2Tj) (i.e. the first term in the right part of equation) is taken zero. 

3 Nonholonomic Diffusion in General Relativity 

In the derivation of relativistic diffusion equations and constructing gravity theories, we have to take into account 
in an appropriate way the fact that the speed of light has a constant maximal value. Geometrically, such a fundamental 
experimental fact is encoded into the special theory of relativity as the condition that the velocity space is a hyperboloid 
(i.e. a special type three dimensional, 3-d, Riemannian manifold) embedded into the 4-d velocity Minkowski space. A 
formal geometric analogy between Euclidian/Riemannian/Finsler etc geometries and respective ones with "pseudo" 
signatures can be preserved by introducing a formal "imaginary" time like in the "early" works on general relativity 

pTnpTi], 

In this paper, we elaborate a geometric and stochastic formalism in order to include in the scheme the relativistic 
diffusion processes with exact solutions for gravitational field equations. Such solutions can be constructed in general 
form only by imposing corresponding nonholonomic constraints on the systems of partial differential equations for 
dynamical and/or stochastic systems. The goal of this section is to consider an extension of the theory of relativistic 
diffusion on (flat) Minkowski and (curved) Einstein spaces when such spacetimes are enabled with conventional h- 
and v-splitting (respectively, with trivial and/or nontrivial N-connection structure). The values of coefficients metrics 
and connections are considered to be given from certain (necessary smooth class) solutions of classical gravitational 
and matter fields equations. In the next section^ the scheme will be completed by elaborating a method of generating 
both classical and stochastic solutions of Einstein equations. 

3.1 The special relativistic nonholonomic diffusion 

The velocity space in special relativity is characterized by a noncompact hyperbolic structure which for the 4-d 
Minkowski spacetime \M is parametrized by a corresponding relation 

- (v 1 ) 2 + (v 2 ) 2 + (v 3 ) 2 + (v 4 ) 2 = -1, (22) 



The same equations can be considered for the probability density (p(r, u) but for the initial condition ip(r = 0, u) = v?( u ) 
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where v a are normalized velocity variables defined for the typical fiber of tangent bundle T( \M). We can intro- 
duce the laboratory time t — tu 1 /c, where c is the light velocity and r is the proper time, and express v 1 — 

1 /2 

[f + (v 2 ) 2 + (v 3 ) 2 + (v 4 ) 2 ] . The local coordinates on T( \M) can be parametrized in the form u a — {x l .y a = v a ), 
where = 1,2,3,4 and a,b,c... = 5,6,7,8. It is considered a conventional 1+3 splitting for \M. We shall 

write h T( \M) if the v-coordinates are subjected to constraints of type (f2"2"j) and say that its typical fiber space is 
a hyperbolic velocity space (when the hyperboloid is embedded into the 4-d Minkovski spacetime). Such a space is 
naturally enabled with a corresponding metric, h^(v c ), and linear connection (Christoffel connection coefficients on 

the hyperboloid), 7~_(v e ), when 

Kb^) = - wV/fr 1 ) 3 , T|(« e ) = v%v (23) 

where (in this subsection) a,b, ... = 2,3, 4. We can prescribe additionally any nonholonomic 2 + 2 and/or 4 + 4 splitting 
on \M and/or T( \M) by prescribing on such spaces corresponding nonholonomic distributions with associated N- 
connection structures of type ((D). A nonholonomic distribution in T( \ M) is given by N = {Ng(u a )}, local coordinates 
u = {u a = (x l , v a )}. The d-metric structure g on T( 3 M)) can be written (in a form similar to ©) 

g=r fe dx i ® dx 3 + h^dv 3 + ivf dx k ) <g> {dv l + N l k dx k ) , (24) 

where rjij = (— , +, +, +), for which, using formulas ([7|). we can compute the corresponding coefficients of the canonical 
d-connection = (L) k M, C%,C$). 

a/3 v J K bk J c be' 

For nonholonomic distributions with Ng(u a ) when J?"^ = 0, see formulas ©, we get C~L = 7|L(w e )- But, in 

general, we work with arbitrary frame and coordinate transforms and the fundamental geometric objects on h T{ \M) 
with nontrivial h- and v-splitting induced by N, when g = { g 3/ g = [gij(u a ), h~^(u a )]}. It is possible to introduce 

N-adapted orthonormalized frame bases Eg' = E 2 ^, (u)ds = E" s , (u)Es for which 

= \9i'i' = Vi>j> , \r v = ^] = E 3 3 ,E^ g s? , (25) 

^/ = ^V^., g ih h n r b , = E%E\, h$, 

when the values with "un-primed" indices are given by the coefficients of d-metric ((24)) and ([23)) . In our approach, 
the frame transform are written with capital letters, E^v, = \E l i , , E°~,], for spaces with pseudo-Euclidian and/or 
hyperbolic fiber/velocity space. The values E % i , parametrize moving frames adapted to the position spacetime \M 
and the values E a ~,(T) state moving velocity frames adapted to the typical fiber of h T( \M). 

Working with "hyperbolic geometric data" g (|M)). r 7 g ~ and E^, ((23)) for a prescribed N-connection structure 

N : TT( \M) = hT( \M) ®vT( \M), we can define a model of nonholonomic relativistic diffusion in special relativity 
(i.e. on total space T( \M)) similarly to the constructions provided in section |2~31 Let us consider the frame bundle 

space F( h T( \M)) with local coordinates f = {w Q ,E% f ,} and identify dx % {r) — v l (r)dT, where r can be interpreted 
as an evolution parameter along the world lines of the particles which can be chosen as the proper time (we can take 
r as the temperature in gravitational thermo-field models etc). The N-adapted relativistic stochastic equations on 
h T( \M) are similar to (fT5|) , 

5u s = E a 3 , (r) o 5W S ' + A a (r)dr (26) 

when the velocity coordinates are subjected to the hyperbolicity conditions ((2"2l and the components of drift d-vector 
are parametrized A s (r) = {A 1 {t) = b l (r) , A u (r) = B u (t)). 

The stochastic differential equations (1261) adapted to the nonholonomic hyperbolic velocity structure split into two 
families (respectively for the position coordinates and for the velocity type coordinates), 

Sx l (r) = E\, (r) o 5W 1 ' + b' (r)dr (27) 
dE i i ,(r) = -L%(r)E k i ,(r)odx k , 

and 

8v z = E\ (t) o <5W a ' + B*(r)dT (28) 
6E%,(r) = -&JT)E%,(r)o6v $ . 
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where a nontrivial B a {r) = F a /rriQ is defined by spacial components of a 4-force F a acting on particles of rest mass 
too. For trivial N-connection structure, the equations ([27]) and ((28)) transform respectively into relativistic stochastic 
equations (2) and (3) proposed in Ref. 8 . In our approach, we can introduce a nonholonomic dynamics in the velocity 
space of special relativity, via corresponding N-connection structure and generalized d-connections. This way we can 
model various theories with restricted/brocken Lorentz symmetry etc. 

For simplicity, we omit here the considerations when from (|2"tJ)) a relativistic theory of (A, L)-diffusion is derived. 

E3 



3.2 Nonholonomic diffusion and gravitational interactions 

There are two classes of relativistic diffusion theories for gravity: 

1. The first one is for modelling relativistic stochastic processes on fixed classical curved spacetimes. We have to 
consider, instead of the flat Minkovski spacetime \M and metric g (|M|) . a (pseudo) Riemannian spacetime V and 
a solution g^ in general relativity (GR). The four-velocity of massive particles in GR satisfies the condition 

g^KK^ = -i. (29) 

(here fi,u = 1,2,3,4), which is a generalization of (1221) (we suppose that in any point u 6 V the hyperbolicity 
condition holds in the typical fiber of TV). By using the orthonormalized moving frames of the pseudo-Riemannian 
manifold we get the same formulas as for relativistic diffusion in special relativity, see details in Section 3 of [5]. 
If the Levi-Civita connection V is changed into the canonical d-connection D (or, for other models, D = C D, 
for instance, the Cartan d-connection) , we get relativistic models on nonholonomic (pseudo) Riemann and/or 
Lagrange-Finsler spaces rnill2j . 

2. The second class of theories is that when d- metrics g = {g^} © are exact solutions of Einstein equations in 
a gravity theory for D, or its restriction to V (with Z J a p = in formulas ([9])), when some coefficients g^j, take 
stochastic values, which results in a relativistic gravitational diffusion of metrics in general relativity (for various 
generalizations). We work with nonholonomic gravitational configurations with associated N-connection splitting 
because in such cases we are able to solve the Einstein equations in general forms and analyze mutual diffusions of 
metrics (the splitting of equations is not possible for not N-adapted constructions with the Levi-Civita connection, 
see section [4]). 

In this section we shall elaborate a model of relativistic diffusion for a nonholonomic (pseudo) Riemannian spaces 
V with prescribed N-connection structure N : TV =/iV®uV, when dimV = 4. The constructions are performed 
for the canonical d-conncction D and d-metrics being solution of the nonholonomic Einstein equations with limits of 
type D| Z ^ -> V. 



3.2.1 N-adapted frames and stochastic equations in GR 

Orthonormalized N-adapted frames on V, e a < = e— a , (u)9„ = e a a , (u)e a , with e a being of type ([5]), can be defined 
by any nondegernerated matrix fields e~ a , (u) and e a a , (u) (and/or, respectively, their inverses, e a a (u) and e™ (u)) 
subjected to the conditions 

Va'P' = S a pe a a ,(u) e^,(u), (30) 

where r\ a 'fi' = (1, 1, — 1, 1), and the space-time orientation of coordinates is chosen in a form which will be convenient 
for constructing exact generic off-diagonal solutions of Einstein equations (see next section). 

We can compute the orthonormalized components v a of a 4-vector v a , v a = e Q Q (u)v a , when rj a ipiv a v@ = — 1 
which can be proven using algebraic relations (|29p and (|3T)]) . By definition, D is metric compatible and we can impose 
the condition D,ge Q Q , = and D^e^, = 0. We get such a transformation law for the canonical d-connection coefficients, 

fV(u) = e/(u) (fV„,(u)+e^e«,(u)) , 

14 Corresponding formulas are similar to (|18p - (|2ip with that difference that the "hyperbolic" small Greek and velocity indices, 
for this section, are with "hats" emphasizing the fact that the nonholonomic diffusion evolution is adapted to the condition of 
constant speed of light. 
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which can be decomposed into h- and v-components using spliting of type e Q Q = (^e^ , e a a J ,&p = (e^efe) and 
T 7 a p = (^jk^tk^jc^bc)- The values r a ^, are analogs of spin connection coefficients E a ^ v , in general relativity; 

pa' v pa' 

L 0v'\Z->O ~> 1 0u'~ 

The N-adapted parallel transport of a d-vector v a with respect to an orthonormalized N-frame epi is defined by 
line elements 

Sv a ' = -f^i/V'fo/, where 5u p = e^,(u a )/'k 

We can say that the change of spacial components (labeled by small Greeck indices with hats, for a 3-d velocity space 
which in the relativistic case is subjected to the condition of hyperbolicity (|22p) of the velocity vector d-field v a is 
driven by a formal (gravitational) force T a ^,{vP). One could be additional contributions from a "stochastic force" 
(such a noise can from any classical or quantum gravitational, or matter fields, fluctuations) associated to a Wiener 
process 8W a . For relativistic constructions, we have to consider stochastic processes along the orthonormal frames 
E~ (v@) as we considered in (1251) . In this subsection, the position space is not the Minkovski spacetime \M but a 
nonholonomic (pseudo) Riemannian one V; we can generalize the constructions considering that in any point u G V 
it is defined a hyperbolic velocity space with the "fiber" metric and connections determined by 

M w7 ) = 5 ap ~ « a ^/(« 3 ) 2 , 7^K) = 

as in (|23[) but with that difference that in this subsection indices run different values (for instance, a = 1,2,4) than 
those in used in subsection 13.11 

Considering that on a nonholonomic (pseudo) Riemannian manifold the Wiener process is moved along the or- 
thonormalized frames E~ (v@) in the 3-d hyperbolic velocity space on u@ £ V, when (with summation on repeating 
indices) 

E%E% = h s \ equivalently h e $ES,E? = 8 &? „ (31) 
where h a P is inverse to the hyperbolic metric h s ^, the infinitesimal motion of the velocity v a is determined by equations 

Sv s = EgSv s ' and 5E%(t) = - 7 L (v?)El,6v$ . (32) 

A random N-adapted curve on V, parametrized in the phase-space can be introduced as a Wiener process: 
dv a — > 8W a . We can consider a formal noise force nmse B a = E%,(t) o 5W a . Following constraints (|3"0)) and 
algebraic relations pip , we conclude that the relativistic diffusion on spacetime V is defined by a more restricted 
system of coefficients because there are admitted only such E~, (r) when a direct relation to the hyperbolic geometry 
is established. We can say that we model via T a a v , a consistent description of Markovian diffusion in general relativity 
but keeping certain " geometric memory" on nonholonomic h-v-splitting and spacetime (pseudo) Riemannian geometry, 
all encoded into a corresponding Stratonovich relativistic calculus. 

The stochastic differential equations describing the N-adapted relativistic diffusion of gravitational and external 
force field^Ifl are written 

8u a = e a a ,(u p )v a 'ST, (33) 
6v s = E%(t) o SW s ' - f% 1 ,{u' 1 )e p a ,{v?)v' 1 'v a ' 5t + ex B s Sr 

where a possible additional external force ex B a = F a /mo is defined by spacial components of a 4-force F a acting 
on particles of rest mass mo and r is a parameter defined along of world line of particles (in our case, moving with 
nonholonomic constraints on V); we can consider r as the phase-space proper time. 

The system of stochastic equations (|32p and (j3"3")l is a respective analogous of (|2T|) and (|2"51) . There are two substantial 
difference between such systems of equations: the first one is for general relativity when the N-connection structure 
is defined for the base spacetime manifold but the second one is for special relativity nonholonomically extended as 
a stochastic geometric model for N : TT( \M) = hT( \M) © vT( \M). Finally, we note that, for instance, there 
are satisfied sufficient and necessary conditions for the existence and uniqueness of N-adapted relativistic stochastic 
differential equations (I33p if the drift and diffusion coefficients are subjected to uniform Lipschiz conditions (see details 
in [23"1I24"II!?5] ) and the stochastic process X(r) = {u(r),v(r)} is N-adapted to the Wiener process W Q (r), when the 
output X( 2 r) is a function of W Q ( 1 r) up to that time, for V < 2 r. 



i.e. the nonholonomic Langevin equations in general relativity 
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3.2.2 Diffusion and N-adapted stochastic GR processes 

A model of general relativistic of (A, L)-diffusion on the fiber bundle F(V) with local coordinates = {u a = 
(x\ y a ), v 13 , E~, } can be derived from (|3"2")1 and (f3"3")l . The N-adapted diffusion operator F ( V )A can be construction 
similarly to (|17p by using operators 



p 



where the force -B e consists respectively from the gravitational and external components, 

B e = -f%,(<)e^,(^V>' + ea; B f . 

For real applications, it is convenient to use the operator P A as the projection of F ( V )A on the phase space 
with coordinates r = {u a = (x l , y a ), v 13 } when for corresponding functions f and / the condition F ^ v ^Af( E r) — 
p A/(u, v). The N-adapted diffusion operator in the phase space is given by 

P A = V A + e a a ,( U P)v a 'e a + B s -^, 
where the Laplace-Beltrami operator in the hyperbolic velocity space is 

"A = 5 & P'E%—^EP —U 
a dv a 0'd v P 

h sp ( & ~ _J d_ ( r— ^ _d_ 

This operator is self-adjoint, "A = ( "A) + . 

The corresponding backward Kolmogorov equation for the N-adapted general relativistic stochastic processes ((33]) 
is written in the form 

— <p(r, u,v)= p A^(r, u,v). 

Using the adjoint of d-operator P A, it is possible to construct the corresponding Fokker-Planck equation in phase 
space (this equation in general relativity is also called the Kramer equation). Introducing the probability density 
function <P :— ip(r, u, v) (as the transition probability <£(u, v, r|uo, Vrj, r = 0)), we write the Fokker-Planck N-adapted 
equation in general relativity on a nonholonomic spacetime with d-metric g a p © and 3-d hyperbolic metric h~^: 

l~a\ 

where the first two terms in the right side are the divergence d-operators in, respectively, the position and velocity 
spaces. 

We use the phase-space proper time r in the N-adapted general diffusion equation ([3"3"f (see detailed explanations 
in |24j how nontrivial torsion terms can be included in the drift coefficients with additional terms in B a ). Alterna- 
tively, we can consider parametrizations in terms of coordinate time (in this work, u 3 = t) which is convenient for 
introducing gravitational and external forces fields. In such cases, the observer time with N-adapted infinitesimal 
element Su 3 = 8y 3 = e 3 a ,(u l3 )v a St is a function of the proper time r and the space and velocity variables. A simple 
analogy with non-relativistic diffusion formula is possible for such frames of references when e 3 3 , = which can 
be introduced in N-adapted Arnovitt-Deser-Misner (ADM). We can introduce 5u 3 = 5y 3 = \/\hzi^u^\v 3 8T, where 

1/2 

v 3 = [l + (v 1 ) 2 + (v 2 ) 2 + (f 4 ) 2 ] . Such a model of general relativistic diffusion in ADM variables, for trivial N- 
splitting is elaborated in [8] (see formulas (36)-(46) in that work). To generate exact stochastic solutions of Einstein 
equations we have to consider nontrivial N-connection structures which makes the stochastic/diffusion theory more 
complex but allows us to separate and integrate the fundamental field equations. 
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4 Exact Stochastic Solutions in Gravity 

The anholonomic deformation/frame method of constructing exact solutions in gravity jT71ll8U2f| can be extended 
to a formalism with stochastic processes and diffusion. In this section, we provide necessary geometric preliminaries 
on gravitational field equations on nonholonomic (pseudo) Riemannian manifolds, show how such equations can be 
formally integrated in very general forms and provide some general criteria/conditions when certain components of 
metrics and connections are induced by stochastic generating functions and corresponding diffusion processes. 



4.1 The Einstein equations on nonholonomic manifolds 

In standard form, the Einstein equations on a nonholonomic (pseudo) Riemannian manifold (spacetime) V are 
written in terms of the Ricci tensor, R p$, and scalar curvature, R, for the Levi-Civita connection V, for a given 
source, i.e. energy-momentum tensor for matter, T^F^l 

R fis - 7^gi3sR = xTps, (35) 

where >t = const. It is not possible to integrate analytically, in general form, this system of partial differential equations 
because of it generic nonlinearity and complexity. In the above mentioned works (see also references therein), we 
proved in details that very general integral varieties can be constructed if we rewrite the equations (|35[) in terms of, 
for instance, the canonical d-connection D, 

R 05 - ~g/w S R = r fiS , (36) 
L^=e„(i\7), C} b =0, 0^ = 0. (37) 

In the above formulas, R p$ is the Ricci tensor for r 7 Q(3 , S R — g^R pg an d T ' ps is constructed for the same metric 
but with D, similarly to formulas with V, when Yps —> xTps for D V. If the constraints (pTT)) are satisfied the 
tensors T 7 Q( g © and Z 1 a ^ from Q are zero and are equivalent to (f55"j). Using nonholonomic deformations and 
"non-tensor" transformation laws for the coefficients of the linear and d-connections, we can satisfy the condition 
r 7 Q/3 = r 7 a/3 , with respect to N-adapted frames §5§ and ([5]), see ©, even V. 

Any exact solution of Einstein equations can be parametrized as a metric (J3J . In classical gravity the coefficients 
gij = diag[g t = ^ , h ab = diag[h a = r) a °h a ] and iVf = w ( = r)f °Wi, N% = n, = nf °n, are certain smooth 
and/or singular "non-random" classes of real functions defining, for instance, a black hole/worm hole /cosmological 
etc solution. The main goal of this paper is to prove that we can construct new classes of Einstein equations, generated 
by some stochastic/random gravitational 'polarizations' rj a and rjf when the nonholonomic deformation of metric 
°g=[ °gi, °h a , °N%] — > ?? g=[ gi,h a ,N%] result both in solutions of equations ((361) and (j37|) . or (|35|) . and any 
type of stochastic/diffusion equation (for various purposes, we can consider certain variants of relativistic diffusion, 
backward Kolmogorov, Fokker-Planck etc equations). For explicit constructions, we can consider that °g is a (pseudo) 
Riemannian metric (which can be, or not, a solution of classical Einstein equations) but impose the condition that 
a metric ''g define a solution (or a class of solutions) of (|56")l when certain coefficients are additionally generated by 
some stochastic/diffusion processes in curved spaces. 

Various types of stochastic generalizations of Einstein manifolds (with mixed types of non-random and random 
processes, mutual diffusion of gravitational and matter fields etc) are called in this work as stochastic Einstein spaces. 
Here we emphasize that because of generic nonlinearity of gravitational field equations the solutions may be with 
chaos, stochastics, fractional behavior etc even we may put certain well defined classical boundary /initial conditions 
on integration functions. Such nonlinear stochastic gravitational and matter field configurations are with a very complex 
and mixed random-sure spacetime structure and very sophisticate rules of stochastic and nonholonomic differentiation 
and integration. 



16 In brief, the gravitational field equations in general relativity are defined geometrically in this form: Denoting by V = {-T^ 7 } 
the Levi-Civita connection (uniquely defined by a given tensor g to be metric compatible, Vg = 0, and with zero torsion), the 
coefficients of necessary tensors are computed with respect to an arbitrary local frame basis e a = (e;,e a ) and its dual basis 
e /3 _ ( e J g 6 ) Using the Riemannian curvature tensor 1Z = {R a p^ e } of V, we define the Ricci tensor, 1Zic = {R ps 4= R a p a s}i 
compute the scalar curvature R = g^ 5 R ps, where g 135 is inverse to g a p. 
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4.2 Generating stochastic solutions of Einstein equations 

For an ansatz of type ([3]), the Einstein equations ([55)1 for D with a general source of typ43 Y a )3 = diag[Y 7 ; Y\ = 

T2 = Y2(x k ,v);Ts = Y4 = Yn(x k )] transform into a system of nonlinear partial differential equations with separation 
of equations for h- and v-components of metric and N-connection coefficients, 

* - « - " ff " % + °" ~ % " - - T ' 1 ^ < 38) 

R3k = 2hl [hi ~ -2hT ~ W + Ahl {— + —)-^hT=°> (40) 
hi »» , ( hi , „ 3 , n£ 



w* = e t In |/i 4 |, e fe Wi = eiWk, n* = 0, ftrifc = dfcn; (42) 

In brief, we wrote the partial derivatives in the form a* = da/dx 1 , a' — da/dx 2 , a* = da/dv. The ansatz and 
resulting system of equations does not depend on variable y 4 (we do not have terms with d/dy 4 , i.e. our original 
ansatz was taken with one Killing symmetry; see references [T7iri8j how to construct exact solutions with general 
"non-Killing" symmetries). The constraints f4"2")l have to be imposed additionally if we wont to satisfy the conditions 
(1371) and generate solutions in general relativity just for the Levi-Civita connection V. 

The above system of equations can be integrated in very general forms |17U18j . in sure variables, by integrating 
step by step the equations beginning with (|38[) . which independent from another ones (relating two coefficients of 
h-metric, g\ and <?2, and source T 4 , all depending on two variables x k ), then (l3l?l) (relating two coefficients of v-metric, 
/13 and hi, and source Y2, all depending on three variables x k ,v, where we take v = t; the time coordinate can be 
related for ansatz of type <j3j) with the parameter r via relation St = yj\h^\8t). The equations (|40l) is an algebraic 
one for N-connection coefficients Wk and we have to integrate two times the equations (|41l) in order to compute the 
N-connection coefficients n^; in both cases, we can generate sure or stochastic solutions, depending on the type of 
solutions we have for /13 and/or /14. Finally, the equations (|42p impose additional constraints on integral varieties for 
([55|) - (f4T|) eliminating induced (by off-diagonal coefficients of metric, i.e. N-connection) torsion and constraining the 
system of equations and solutions to be just for V. 

The h and v separation allows us not only to generalize the constructions for stochastic metrics and N-connections 
but also to analyze the conditions when a non-random h-metric structure model gravitational diffusion processes in 
the v-subspace. 



4.3 Stochastic solutions with h$ 4 ^ and T 2 4 ^ 

We consider a metric ([3]) when 

''g = e^ x " ) dx i ® dx l + h 3 (x k ,t)e 3 ®e 3 + h 4 {x k , t)e 4 (g>e 4 , (43) 
e 3 = dt + w i (x k ,t)dx i ,e 4 = dy 4 + n i (x k ,t)dx i 

is supposed to be a solution of (|3"8"l) - (f4"Tj) with gi — g2 — ) defining a non-random h-metric, as a solution of 2-d 
Laplace equation (the h-components of the Einstein equations transform to this simple equation which can be used 
for generating sure, or random solutions) 

$ + iP" = 2Y i (x k ), (44) 

and /13, defining a stochastic v-metric (in result, wt and/or rij can be also stochastic solutions). 
Introducing values 

<P = ln\ ai = hld i( f>, (3 = hl0*, 7= (ln|/M| 3/2 /l^ 3 |)\ (45) 

17 parametrization of energy-momentum tensors in the above presented form are possible by corresponding nonholonomic 
frame and/or coordinate frame transform for various types of matter sources, including some general and important cases with 
cosmological constants and various models of locally anisotropic fluid/scalar field/ spinor/ gauge fields interactions on curved 
spaces 
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the equations (j4"0")) . (|4"Tj) are respectively written in the form 

pwt + on = 0, (46) 
<*+ 7 < = (47) 

The type of solutions for N-connection coefficients depend explicitly on the type of solutions we construct/choose for 
the v-metric. 



4.3.1 Non-random/sure solutions 

For sure coefficients, the equation transform into 

hl = 2h 3 h 4 r 2 {x i ,t)/<P*. (48) 

If h\ ^ 0;T2 ^ 0, we get <f>* ^ 0. Prescribing any non-constant cj> = (j)(x l ,t) as a generating function, we can construct 
exact solutions of (|48[) -(|47 p . Integrating on t, in order to determine /13, /14 and rii, and solving algebraic equations, for 
Wi, we get 

ft3 . ± !£&2L , ht - ± 2 / <-' 2 f ■ '"'V ,49) 

^2 J ^2 

ttfi = «i = l n k (x*) + 2 n k (**) f[h 3 /(^\h^\) 3 }dt, 



where °/i4(x fe ), 1 n k (a? 1 ) and 2 n k (x l ) are integration functions. We have to fix a corresponding sign ± in order 

to generate a necessary local signature of type (+ H h) for some chosen 0, T2 and T4. Such solutions include as 

particular cases the classes of solutions for a nontrivial cosmological constant X,; = A, or nonholonomic configurations 
with polarizations of such constants, A — > h X(x k ) = Tn(x k ) and A — > v X(x k ,t) — T2(x k ,t). 

In Refs. |171ll8ll2"T] . we studied various classes of parametric and non-parametric exact solutions with coefficients 
of type (|44|) and (|48|) describing black ellipsoids, locally anisotropic wormholes and cosmological solutions and their 
noncommutative generalizations. Those constructions can be generalized for stochastic Einstein spaces if <p and/or 
Ti are taken to be certain stochastic functions. 

4.3.2 Metrics with h-diffusion of generating functions 

This class of solutions if we chose, for instance, a generating function 

4>{x k , t) -> 0{x k , t) = <j)(x k ,t) + m$(x k ,t), (50) 

where (f>(x k ,t) ~ f(r, u) = f(T,x l ) is defined as as diffusion process from h-space to v-space by a corresponding 
backward Kolmogorov equation for the N-adapted general relativistic stochastic processes (|2"0)) and/or (|2"Tj) (in order 
to state the solutions in exact form, we consider only h-operators and fix A' 3 = in V A (|18p). where r t, 
v = vo = const and the Laplace-Beltrami d-operator A (|T5)) is computed as A = hA for g^j = &ye'^ x \ We can 
consider roasa real parameter which for w = transforms our system into a non-random one. For solutions with one 
Killing symmetry such parameters can be always introduced. Nevertheless, there are substantial differences between 
the parameter w and those for sure solutions with superposition of Killing symmetries and/or with noncommutative 
parameters 9, see [H]. In this work, a nontrivial w results in stochastic behavior of (some) coefficients of metrics and 
connections, i.e. of gravitational fields. 

The generalized Kolmogorov backward equation (describing diffusion of gravitational h-components on r, i. e. into 
^-components) is 

cVf (r, x l ) = ftA/((r, x*) = &A/((t, x l ), (51) 
l(0,x l ) = f(x i ). 

The corresponding generalized Fokker-Planck equation is 

d T F = ZhAF, (52) 
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where F = F(t, 1 x l ;0, 2 x l ) is the transition probability with the initial condition F(0, 1 x % ;0, 2 x l ) = S( 1 x % — 2 x % ) 
for any two points 1 x l , 2 x l G V and adequate boundary conditions at infinity. The same equations can be consid- 
ered for the probability density ip(r, x l ) but for the initial condition ip(r = 0,x l ) = °(p(x l ). Stochastic components 
(f>(x k ,t) in (|49|) and ((50|) generated as diffusion solutions of (l5Tj) and/or ((52|) provide some "simple" and explicit 
examples of stochastic metrics in Einstein gravity, all distinguished by corresponding nonholonomic constraints. Un- 
der general frame/coordinate transforms, the parametizatons change into mixed holonomic/nonholonomic and non- 
random/stochastic variables both on h- and v-subspaces. More general types of nonholonomic diffusion equations can 
be also considered but the solutions will be generated in non-explicit forms. 

Fixing 2 7ifc =0 in (HHJ), we generate non-random solutions for the N-connection components n^. But Wi will 
have positively a stochastic character, as well both /13 and /14, if the generating function <fi(x k , t) contains a nontrivial 
diffusion component. Such random gravitational processes can be induced even as vacuum gravitational configurations 
for smooth coefficients of Einstein equations 



4.3.3 Solutions induced by random sources 

There is an alternative possibility to generate random gravitational processes than those with stochastic gener- 
ating functions. This can be seen from formulas (|48[) and (|49l) for certain random/noise behavior of source T^. Such 
constructions are similar to those for stochastic gravity (based on the Einstein-Langeven equation with additional 
sources due to the noise kernel, see a review and references in [53]). Nevertheless, it should be mentioned here that 
our anholonomic deformation method is different from the methods of deriving solutions in stochastic gravity. Via 
nonholonomic transforms/deformations we can generate in a "non-perturbative" manner different classes of exact 
solutions with generic off-diagonal terms and mixed holonomic-nonholonomic variables and matter sources. 

fn a particular case, we can consider that a stochastic T2 can be constructed as a random vacuum polarization 
of gravitational constant, when vacuum gravitational fields are nonholonomically distorted as a diffusion process, for 
instance, modelled by solutions of equations (|5 1 [) and/or ([5^)) associated to some effective matter fiels/cosmological 
constant. 



4.3.4 Diffusion to the Levi-Civita conditions 

In order to construct exact solutions for the Levi-Civita connection, i.e. of standard form of Einstein equations 
written with respect to N-adapted frames, we have to constrain the coefficients (|4T)|) of metric (|43[) to satisfy the 
conditions (|42|) . There are various possibilities. For instance, we may consider a classical sure metric (which may 
be, or not, a solution of gravitational field equations with nontrivial distortion) and constrain the system in such a 
way that gravitational diffusion will result in stochastic vacuum or Einstein spaces. There are scenarios when a sure 
Einstein configuration is stochastically transformed into nonholonomic configuration with nontrivial distortion of of 
the Levi-Civita connection. 

The above mentioned gravitational diffusion evolution models depend on the class of additional constraints we 
impose on generating and integration functions. For instance, we can chose that 2 ri]~ (a: 1 ) = and 1 rik (x l ) are any 
functions satisfying the non-random conditions di 1 rik = c\. 1 rii. The constraints for <p(x k ,t) can be for random, sure 
and/or of mixed nature variables, following from constraints on N-connection coefficients Wi = ~di<j)/(jf , 

(wMY + vhW {hM)* + dM<f>] = 0, 

di w k [cj>] = d k Wi[4>], (53) 

where, for instance, we denoted by h^ifi] the functional dependence on <fi. So, if <fi is stochastic, we have to consider 
that are some relations on mathematical expectations etc. Such conditions are always satisfied for cosmological 
solutions with <j> = <fi(t) or if <j> = const (in the last case Wi(x k ,t) can be any non-random and/or stochastic functions 
as follows from |46|) with zero /3 and a,, see l[4"5"])). 



4.4 Special cases of stochastic Einstein spaces 

We can construct such solutions for certain special parametrizations of coefficients for ansatz (|4T)|) subjected to the 
condition to be solutions of equations (|44l) - (|47[) and certain nonholonomic diffusion equations on curved spaces. 
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4.4.1 Vacuum gravitational diffusion with h\ = 

The equation (jSHJ) can be solved for such a case, h\ = 0, only if T 2 = 0. So, in N-adapted frames, the v-components 
of gravitational equations are vacuum ones. We can consider any functions Wi(x k ,t), being sure or random ones, as 
solutions of P0|) , and its equivalent (|4^|) . because the coefficients (3 and Ofj, see (f4"5"j). are zero. 

We find nontrivial values of by integrating (|47[) for h\ — and any given h 3 which results in ni — 1 ri k (a; 1 ) + 
2 n k (x l ) f h 3 dt. Choosing h 3 to be a stochastic function, we have to compute J h 3 dt as a Stratonovich stochastic 
integral, when f h 3 dt — > J h 3 o dr. We can consider any 51 = 32 = \ with ip(x k ) determined by ([44)) for a given 
T4(x fc ); for simplicity, we can consider only non-random solutions for the h-metric. 

Summarizing the constructions, we get a class of stochastic gravitational solutions defined by ansatz 

"g = e^ {xk) dx l ® dx l + h 3 {x k ,t)e 3 ®e 3 + °h 4 (x k )e i ®e i , (54) 
e 3 = dt + Wi(x k ,t)dx l , 

e 4 = dy A + [ x rifc (a; 4 ) + 2 n fc (a;') J h 3 dt -> / fr 3 o dr]dx fc , 

for arbitrary generating stochastic and/or sure functions h 3 {x k , t), Wi(x k , t), °h 4 (x k ) and integration functions 1 rifc (a; 4 ) 
and 2 n,t (a; 4 ) . 

The conditions (|42|) selecting from (|54|) a subclass of solutions for the Levi-Civita connection transform into the 
equations 

2 n k (V) = and 9, = <9 fc n», 

w* + di hi = and di w k = d k Wi, 

for any such Wi(x k ,t) and fi4 : (x k ). This class of constraints do not involve the generating function h 3 (x k ,t). So, we 
can consider sure constraints to the Levi-Civita configurations for a stochastic h 3 . In general, we can model nontrivial 
Levi-Civita diffusion processes for sure and/or stochastic h 3 , but constraining the gravitational diffusion process to 
conditions h\ — and T 2 = 0. 



4.4.2 Gravitational diffusion with h% = and h% ^ 

Such stochastic spacetimes are defined by ansatz of type 

v s = e^dx 1 ® dx l - °h 3 (x k )e 3 ®e 3 + h 4 (x k ,t)e 4 (E)e i 7 

e 3 = dt + w i (x k ,t)dx i ,e 4: = dy 4 + n t (x k , t)dx\ (55) 

where gi= g% = e^ x> °\ with ij](x k ) being a solution of (|44l) for any given T 4 (x k ). The function h 4 (x k ,t) must satisfy 
the equation flU} which for h 3 = is just 

h?-^-2°h 3 h i r 2 (x k ,t) = o. 

The N-connection coefficients are 

Wi = -d^/ft, rii = l n h {x l ) + 2 n k (x*) J [1/( V / M) 3 ]^, 

when <j> = In \h\\/ */] h 3 li4\\. Fixing °h 3 — 0, we can always eliminate possible stochastic contributions from T 2 with 
a sure value /14 being a solution of h** = (h\) 2 /2h 4 . This impose us to consider classical non-random values for <fi 
and, as consequences, for Wi and rij. 

If we consider any °h 3 / 0, the stochastic gravitational processes will be induced by "noise" in T 2 . So, the stochastic 
metrics of type ([5"5"j) are generically defined by a stochastic source T 2 . Such nonholonomic stochastic gravitational 
configurations are very different by those described by ansatz (|54p when the gravitational diffusion is of generic 
vacuum type. 

The Levi-Civita configurations for (|55|) are selected by the conditions (j42]) which, for this case, are satisfied if 

2 n k (x l ) = and di 1 n k = d k 1 n i , 
and (wi[(j)\j + Wi [4>] (ji 4 [(j)\j + dihi[<f>] = 0, 

di w k [(f>] = d k Wi[<f>]. 
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Such conditions are similar to (|53l) but for a different relation of v-coefficients of metric to another type of generating 
function (f> which can be stochastic only for random values of Ti . They are always satisfied for cosmological solutions 
with <fi = <p(t) or if (f> = const. In the last case Wi{x k ,t) can be any stochastic or sure functions as follows from (|46l) with 
zero p and on, see P5|) . So, for some special configurations, random values of w k can be induced via nonholonomic 
vacuum gravitational diffusion even the source T2 is constrained to be non-random. 

4.4.3 Solutions with constant generating functions 

Fixing 4> — 4>q = const in (|45|) . with ft| ^ and h\ ^ 0, we can express the general solutions of (j44"]) -(l47" l) in the 
form 

n g = e ^ k )dx l ®dx l - °h 2 [f* (x\t)] 2 |sr (x\t) |e 3 (g>e 3 + f (x\t) e 4 ®e 4 , 

e 3 = dt + Wl {x k ,t)dx l , e 4 = dy 4 + n k (x\t) dx\ (56) 

where °h = const, g\ = g% = e^ x \ with iji{x k ) being a solution of (f4~4"| for any given T4(x k ), and 

,r (x\t) = Q[0] -M/ r 2 (x k ,t)[f (x\t)] 2 dt. 

In a stochastic sense, we should consider J ...dt J Stratonovich St. 
The N-connection coefficients iV 3 = Wi(x k ,t) Nf — ni(x k ,t) are 

djtfr jx k ,t) 

n k = V (x*) + 2 n k [ ^.^J & (x*,t) dt, 

J \f{x\t)] 2 

with necessary generalizations to Stratonovich stochastic integrals. 

We must take <?4[o] (x 4 ) = ±1 if qr (x % ,t\ = ±1 for Y2 — > 0. In such a case, the functions /13 = — °h 2 [/* (x l ,t)~\ 

and hi = f 2 (x l ,t) satisfy the equation (|48p written in the form ^/l^sl = Q h{y/\h^)*, which is compatible with the 
condition 4> = 4>q. 

The gravitational diffusion for this class of solutions has two sources: the first one can be if / (x l , f J is a random 
generating function and the second one is for random sources T2(x k , f). Gravitational diffusion mix such contributions 
even we distinguish them nonholonomically for certain configurations. 

The subclass of solutions for the Levi-Civita connection with ansatz of type ([55| is selected via conditions (l4"2l) . 
We can chose that 2 n k (x l ) — and 1 n k are any functions satisfying the conditions dt 1 n k — d k 1 ni. The 
constraints on values Wi = —di^r/^j- result in constraints on $r, which is determined by T2 and /, 

(wifcr])* + Wi[sr] (h^r})* + ft/i 4 [?r] = 0, 

di w k [<;r} = d k w^r], (59) 

where, for instance, we denoted by /i4[<rr] the functional dependence on <;x including random and non-random contri- 
butions both from vacuum gravitational configurations and stochastic processes for matter. Such conditions are always 
satisfied for cosmological solutions with / = /(f). For D, if Y2 = and cf> — const, the coefficients Wi(x k ,t) can be 
arbitrary functions (we can 6x^ = 1, which does not impose a functional dependence of Wi on qx) & s follows from (|46[) 
with zero j3 and on, see (|45[) . As in the previous case, such N-connection components can be stochastic ones as arising 
from some vacuum gravitational configurations. To generate solutions for V such n>i must be additionally constrained 
following formulas (|59|) re-written for «Jj[5r] Wi{x k ,t) and /14 [<rr] h 4 (x l ,t) . Such scenarios are typical ones with 
nonholonomic gravitational diffusion from/to Levi-Civita configurations. 



(57) 
(58) 



4.4.4 Non-Killing stochastic solutions 

We note that any sure and/or stochastic solution g = {g a 'P'{u a )} of the Einstein equations (l3"6")l and/or (1331) with 
Killing symmetry d/dy (for local coordinates in the form y 3 — t and y 4 = y) can be parametrized in a form derived in 

this section. Using frame transforms of type e a = e a a e a ', with g Q( g = e a a e^pg a 'p', for any g Q( g we relate the class 



Sergiu I. Vacaru: Nonholonomic Relativistic Diffusion and Exact Solutions for Stochastic Einstein Spaces 



19 



of such (inhomogeneous) cosmological solutions, for instance, to the family of metrics of type ((321) EE in explicit form, 
the solutions can be modelled as explicit diffusions from the h- to v-componcnts of geometrical/physical objects. 

Following our recent results on constructing general solutions in Einstein gravity and modifications |171I181[2"T] . 
we can construct 'non-Killing' solutions depending on all coordinates. Such general classes of sure solutions can be 
parametrized in the form 

g = +g i (x k )dx i ®dx i + cu 2 (x 3 ,t,y)h a (x k ,t)e a ®e a , 

e 3 = dy 3 + Wl (x k , t)dx\ e 4 = dy A + n t (x k ,t)dx\ (60) 

for any w for which 

e k uj = <9fcW + w k uj* + n k du}/dy — 0, (61) 

when ((BDJ) with u? = 1 is of type (J3J). 

A class of stochastic solutions with nontrivial sure factor cj satisfying the sure conditions (I5T1) with sure values 
w k and rik but stochastic h a can be generated similarly to the case of nonholonomic sure Killing configurations. In a 
more general approach, we can impose the conditions (|61[) for any stochastic values uj^n k ,Wk and h a . With respect 
to N-adapted frames such conditions can be satisfied for certain particular, explicit examples, configurations or in 
non-explicit forms. 



5 Summary and Conclusions 

In this work, we have studied the best strategy, and elaborated a geometric method, for generating stochastic exact 
solutions of Einstein equations in general relativity and modifications. 

Our proposal is to follow respectively certain key steps and sub-steps for constructing deterministic and stochastic 
solutions for nonholonomic gravitational and matter fields interactions: 

1. Nonholonomic splitting & formal integration of Einstein eqs 

— Let us introduce a nonholonomic 2 + 2 distribution stating a nonlinear connection (N-connection) structure 
N : TV =hV © v~V with associated N-adapted frames, e a = (ei,d a ) and eP = (e 4 ,e a ), on a (pseudo) 
Riemannian manifold V enabled with metric structure g, when d a — d/dy a and e l = dx l for local coordinates 
u fi = (x j ,y a ), or u = (x,y). 

— We adapt all geometric constructions to N and define the canonical distinguished connection(d-connection) 
D = V + Z, where V and Z are respectively the Levi-Civita connection V and the canonical distortion 
distinguished tensor (d-tensor), all defined in a unique metric compatible form by g. Any N-connection structure 
N = {iV?(u)} states a conventional horizontal-vertical (h-v) decompositions of geometric object, for instance, 

D = hD © vD and g =hg © vg, for hg = {gij} and vg = {h a f,} with respect to N-elongated frames, when 
e 4 = d l - iV| d b and e a = dy a + ATf dx l . 

— We write the Einstein equations (|36p and ([57)1 for the geometric data (g, N,D); the gravitational filed equations 
are equivalent to those in standard variables (g, V), see (|35p . for correspondingly defined energy-momentum 
and distortion sources. 

— For a general ansatz, we get that the Einstein equations are equivalent to a system of nonlinear PDE (|3"5 )) - pTj) 
and (|42l) , with a " consequent splitting" of equations which allows us to construct general off-diagonal solutions 
in sure variables; in local coordinate frames such metrics are parametrized in the forrrF^I 

gi + o; 2 (w 1 2 /i3 + w 2 (n 1 2 /i4) ui 2 (w 1 w 2 h^ + n 1 n 2 h4) uJ 2 1x^/13 w 2 n-Ji/± 
ui 2 (w 1 w 2 h3 + n 1 n 2 /i4) g 2 + uj 2 (w 2 2 h3 + n 2 2 h/C) 0S 2 w 2 h^ uJ 2 n 2 Ji4 
lu 2 w 2 w 2 hz /13 

uj 2 n x h^ uj 2 n 2 hi hi 

where, in this work, y a — (y 3 — t,y 4 = y) and spacetime signature is chosen (+,+,—,+). The coefficients 
gk(x z ), h a (x\ t), Wk(x l , i), nk(x l , t), q(x\ t) and ui(x\t,y) can be defined in explicit form by integrating and/or 
differentiating some generating functions. For instance, the class of certain solutions §3§ is with functional 

18 We have to solve certain systems of quadratic algebraic equations and define some e a a {u^), choosing a convenient system 
of coordinates u a — u a (it' 3 ). 

19 as a matter of principle, any solution of gravitational field equations with certain general matter fields sources can be 
represented in such a form by corresponding frame and coordinate transforms; we have to involve certain additional physical 
considerations, suppositions on symmetry of interactions and boundary conditions in order to model realistic gravitational 
interactions etc. 
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dependence on a generating function <f>{x ,£), when h a = h a [x l ,<p] and (following from relations for v-metric) 
Wi — Wi[x l , t, </>] and i%i = rii\x % ,t, <j>\. The type of chosen generating and integration functions (smooth class, 
sure or stochastic character), prescribed symmetries and topology of interactions, boundary/limit conditions 
etc distinguish in explicit form a geometric/physically important class of solutions. 

2. Elaborate a nonholonomic stochastic calculus and related general relativistic diffusion theory on V 

— The theory of stochastic processes and diffusion on curved manifolds, Riemann and Riemann-Cartan spaces 
E3> and bundle spaces is a well developed direction in modern mathematics P^l 

— In this work, we have shown how the mathematical formalism for the relativistic stochastic theory and diffusion 
can be adapted for nonholonomic manifolds with nontrivial N-connection structure. In brief, the constructions 
are those for (pseudo) Riemannian manifolds but with corresponding generalizations with respect necessary 
classes of orthonormalized N-adapted frames, modified linear connections and Laplace-Beltramy operators. 

— The geometric data (g, N,D) and/or sure general solutions are considered as a fixed nonholonomic background 
for " non-integrable" rolling of stochastic processes in tangent bundles with h- and v-splitting. So, at this stage 
the generating/integration functions for the solutions of Einstein equations are taken to be sure functions even 
fluctuations to nonholonomic Einstein-Langeven systems, induced by quasi-classical quantum fluctuations of 
matter, can be considered in a self-consistent form as in so-called stochastic gravity [26]. 

3. Consider random generating functions and sources for formal solutions of Einstein equations 

— The key idea of this work is to consider (pseudo) Riemannian metrics with coefficients generalized in a form to 
include random variables induced by generating functions (for some particular cases, we can consider directly 
some N-adapted coefficients of v-metric, h ab , and N-connection, N?) defined, for instance, by some gravitational 
diffusion processes. 

— Transitions from a "sure" nonholonomic (pseudo) Riemannian configuration ( °g, °N, °D|z=o — ^ °V) to a 
stochastic nonholonomic Einstein spacetime ( v g, V 'N, r, D| 2 ' = o — > V V) are modeled by metrics of type ([3]), 
when gij = diag[gi = r/i °gi\ and h ab = diag[h a = r\ a °h a ] and N% = w ( = r)f °Wi and N% = m = rjf °n, 
are constructed for some sure and/or stochastic gravitational ^-polarizations. The polarizations are defined in 
such a way that ''g generates a class of exact solutions of the Einstein equations with r\ a and 77" being related 
to some stochastic/diffusion processes (as solutions of some generalized and/or nonholonomically constrained 
Kolmogorov/ Fokker-Planck equations) rolled on a "sure" background ( °g, °N, °D). We have to change the 
usual "sure" integrations into stochastic Stratonovich ones (which is the most convenient for curved spaces), with 
possible re-definition of results for Ito stochastic integrals. For real physical applications, we have to introduce 
orthonormalized N-adapted frames related to some certain degrees defined by °g and sure deformations to 
''g but consider that additional deformations are driven by some Wiener processes <5W° (for physicists, all 
necessary concepts on stochastic calculus and diffusion are summarized in Ref. [2"2"]V 

— In explicit form, the nonholonomic transitions between sure and stochastic gravitational configurations are 
modelled by changing of a "sure" generating function into a random one, cf>(x k ,t) — s> <fi(x k ,t) — <f>(x • , t) + 

vatfi(x k ,t), see ptj]) . Such mutual evolutions of deterministic and random phases are possible even for "sure" 
boundary/initial conditions, in classical gravity because of generic nonlinear character of Einstein equations. 
The stochastic behavior may arise for very small values of parameter w and vanish if w — 0. Geometrically, 
we can impose such nonholonomic splitting into conventional h- and v-subsystems with distinguished sure and 
stochastic variables, which allows us to split the equations and generate exact solutions. 

We emphasize here that it would not be possible to elaborate such, in general, non-perturbative methods of 
constructing exact solutions in gravity with "sure" and/or "stochastic" variables, and their mixture, if we do not apply 
corresponding methods from the geometry of nonholonomic manifolds/bundles and N-connection formalism (originally 
proposed in Finsler and Lagrange geometry and further re-considered in Einstein gravity and modifications). 

Let us compare our methods to an alternative approach in stochastic gravity [26] . For various applications in gravity 
and astrophysics, there are considered semiclassical generalizations of Einstein equation with sources (T a b[g])' ren 
computed for suitably renormalized expectation value of the stress tensor operator. Such a renormalized value is 
computed with the scalar field operator satisfying the Klein-Gordon equation on the perturbed metric g a b + hab and 
stochastic source We get the so-called Einstein-Langevin equation, 

1 G ab [g + h] = k{ 1 T ab [g + h]) ren + x£ ab [g], 



the last ones are with nontrivial torsion; torsion can be also induced nonholonomically on Minkowski and Riemann spaces 
but there are canonical anholonomic transforms to the Levi-Civita connection with zero torsion 

21 During last two decades the approach was extended for (in general, suppersymmetric) Lagrange-Finsler spaces, and higher 
order generalizations, and on (pseudo) Riemannian/Lorentz manifolds with various applications of relativistic diffusions in 
modern cosmology and astrophysics. 
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where the left superindex means that only terms linear in the metric perturbations are kept (the stochastic source 
£ a b is regarded to be of the same order as h ab ). The Gaussian stochastic source £, ab is completely determined by the 
following correlation functions: 

(£ab[g;x))t = 0, (€ab\9;x)€ ab \g;y))t = N abcd (x 7 y) = )^{{t ab [g\x),t cd \g\y)}), 

where (. . .)^ denotes the expectation value with respect to the stochastic classical (i.e. non-quantum) source £ 6[ff]> The 
operator t ab [g\ is defined as t ab [g] := T ab [g\ — {T ab [g\) and the bitensor N abc d(x,y), which determines the correlation 
function of the stochastic source, is computed using the scalar field operator satisfying the Klein-Gordon equation for 
the background metric g ab . The bitensor N abc d{x,y) is called the noise kernel, it describes the quantum fluctuations 
of the stress tensor operator and is positive-semidefinite. 

Comparing our approach and the stochastic gravity scheme we conclude that we can include that scheme into the 
anholonomic deformation method considering a subclass of perturbative "fluctuations" of metrics defined by "noise" 
in energy-momentum tensors. In general, following stochastic nonholonomic geometric constructions we can model 
generic nonlinear stochastic processes both in vacuum and non-vacuum gravity. 
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